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Abstract 

In the context of the Type I seesaw mechanism, we carry out a systematic study 
of the constraints that result from zeros in both the Dirac and right-handed Majorana 
neutrino mass matrices. We find that most constraints can be expressed in the stan¬ 
dard form with one or two element/cofactor zeros alone, while there are 9 classes of 
nonstandard constraints. We show that all the constraints are stable under one-loop 
renormalization group running from the lightest right-handed neutrino mass scale to 
the electroweak scale. We study the predictions of the nonstandard constraints for the 
lightest neutrino mass, Dirac CP phase and neutrinoless double beta decay. 
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1 Introduction 


Many Standard Model (SM) extensions predict the existence of right-handed (RH) Ma- 
jorana neutrinos, which are often used to explain neutrino masses via the Type I seesaw 
mechanism [1]. The effective mass matrix of the light Majorana neutrinos is given by 

M = -MdMr-^MI , (1) 

where and M/j are the mass matrices for the Dirac and RH Majorana neutrinos, respec¬ 
tively. 

In the basis where the charged lepton mass matrix is diagonal, the light neutrino mass 
matrix is described by 9 real parameters and can be written as [2] 

M = I/*diag(mi, 7712, , (2) 

where the rrii are real and nonnegative, V = 17diag(l, and 
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After decades of neutrino oscillation experiments, hve parameters (three mixing angles 
and two mass-squared differences) in the neutrino sector have been determined to rather 
good precision. However, the hierarchy of the three light neutrino masses is still unknown, 
and there are two possibilities: fij mi < m 2 < m 3 for the normal hierarchy (NH) and fiij 
m 3 < mi < m 2 for the inverted hierarchy (IH). The remaining four unknown parameters 
may be taken as the lightest mass {mi for the NH, or m 3 for the IH), the Dirac CP phase 
S and two Majorana phases (0i and 02)- The Dirac phase and the mass hierarchy will be 
measured in future neutrino oscillation experiments, and the lightest mass can in principle 
be determined from beta decay experiments and cosmological observations. 

Of the many neutrino mixing models proposed in the literature, texture zero models [3] 
are attractive because of the simple relations on the matrix elements that are found in these 
models. Texture zero constraints can be imposed on the light neutrino mass matrix directly. 
However, in the context of the Type I seesaw mechanism, since the light mass matrix is 
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a product of the Dirac mass matrix and the RH neutrino mass matrix, it is natural to 
consider texture zeros in the Dirac and RH neutrino mass matrices [4], Zeros in the Dirac 
and RH neutrino mass matrix can be realized from discrete Zat symmetries with suitable 
scalar singlets [5]. 

In this paper, we carry out a systematic study of the constraints that result from various 
zero textures in both the Dirac and RH Majorana mass matrices in the context of the 
Type I seesaw mechanism. We hnd that most cases lead to the standard constraints which 
can be expressed in the form of just one or two element/cofactor zeros, which have been 
studied extensively in the literature [6, 7, 8, 9, 10, 11, 12], However, there are also 9 classes 
of nonstandard constraints, i.e., constraints which are not described by simple texture or 
cofactor zeros in the light neutrino mass matrix. We first show that both the standard and 
nonstandard constraints are stable under one-loop renormalization group equation (RGE) 
running from the lightest RH neutrino mass scale Mi to the electroweak scale Mz- Then we 
study the phenomenological implications of the nonstandard constraints using recent data 
measured by neutrino oscillation experiments. We hnd that some cases are excluded, and for 
the rest we obtain preferred values for the lightest mass and Dirac CP phase, which will be 
probed in the next generation of neutrino experiments. We also study neutrinoless double 
beta decay for the nonstandard constraints. 

This paper is organized as follows. In Section 2, we explore the constraints of zero 
textures in the seesaw mechanism. In Section 3, we discuss RGE effects on the constraints. 
In Section 4, we study the phenomenological implications of the nonstandard constraints 
and summarize our results in Section 5. 


2 Classes of constraints 

We consider three generations of RH Majorana neutrinos and we assume Mr is not singular. 
Since the Majorana mass matrix is symmetric, there are 6 independent matrix elements in 
Mr, and 9 independent elements in Mr,. For a case with Nr zeros in Mr and Nr, zeros in 
Mr, there are shonld be at least N — 6 complex constraints on the elements of M, where 
N = Nr + Nr is the total nnmber of zeros element. This may be nnderstood as follows. 
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We start with 15 = 6 + 9 nonzero complex parameters in and respectively. There 
are a total of N zeros in Mr and Mr, so that leaves 15 — iV independent complex parameters. 
We can always insert a diagonal matrix between Mr and Mr (and the same matrix between 
Mr and M^), which effectively absorbs 3 complex parameters (or, putting it another way, 
there are 3 relative normalizations that are redundant in the prediction of the light neutrino 
mass matrix from the seesaw mechanism). This leaves us with 12 - N independent complex 
parameters. In M there are 6 independent elements, so the net number of constraints on 
the elements of M is 6 — (12 — N) = N — 6. In terms of real parameters, we have 24 — 2N, 
from which we can subtract 3 unphysical phases, leaving us with 21 — 2N physical real 
parameters in M. Since there are 9 real parameters in the neutrino sector and 5 of them are 
well determined, in general, if iV > 9, the model will be over constrained and if iV < 6, the 
model will be unconstrained. 

This counting rule only gives the minimum number of constraints. There are accidental 
constraints that result from some special combinations of texture zeros in Mr and Mr, e.g., 
if Mr has a column or a submatrix composed of all zeros, then the constraint det M = 0 
always exists regardless of the counting rule. So in order to obtain all the constraints for a 
given zero structure of Mr and Mr, we have to make sure all of the independent constraints 
are obtained. Our analysis proceeds as follows. 

We scan all possible texture zero structures for Mr. Since M is invariant under a permu¬ 
tation of the columns and rows of Mr with a simultaneous permutation of the corresponding 
columns of Mr, we only consider one representative structure for each set of Mr that are 
equivalent under permutation. For each combination of Mr and Mr, we check all the el¬ 
ements and cofactors, and the determinant of M using Eq. (1) and look for independent 
constraints. We write the constraints in terms of the elements, cofactors or determinant of 
M in the simplest form. The experimental data does not allow more than two zeros in the 
elements or cofactors of the light mass matrix, hence we only consider cases that give less 
than three such zeros. The cases that give a block diagonal mass matrix are also ignored 
because they lead to a zero mixing angle, which is excluded by current experimental data. 
For the other cases, we analyze the elements of M, and make sure we obtain all of the inde¬ 
pendent constraints. We hnd that most constraints can be expressed in the standard form 
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Table 1: Classes of the nonstandard constraints that can not be expressed in terms of one 
or two element (Mq,^) or cofactor (Ca/s) zeros alone. Here a /3 ^ ■y, and no sum is implied 
in the constraint equations. 


Class 

Nr + Nr, 

Nonstandard constraints 

No. of constraints 

1 

4 + 5 

Maa = Mpp = 0 and det M = 0 

3 

2 

4 + 4, 3 + 5, 3 + 4, 2 + 5 

Meta = 0 and det M = 0 

3 

3 

3 + 5 

Map = 0 and det M = 0 

3 

4 

4 + 4 

Maa = 0 and MppMa^ = 2MapMp^ 

6 

5 

3 + 5 

Ca.oL 0 Rlld C 

6 

6 

4 + 3, 3 + 4, 2 + 5 

MaaCaa = det M 

3 

7 

4 + 3 

MaaCaa = ^MapMa^Cp^ 

3 

8 

3 + 4 

Mpp det M = -MlpCaa 

6 

9 

4 + 3, 3 + 4, 2 + 5 

det M = 0 

1 


that can be written in terms of just one or two element/cofactor zeros, but there are also 
nine nonstandard constraints, which are listed in Table 1. The structures that lead to the 
nonstandard constraints are provided in Appendix A. Classes 2, 3, 4 and 5 are discussed in 
Ref. [4], with constraint equations in a different (albeit equivalent) form. Classes 1 and 9 are 
subsets of well-known models with the additional constraint det M = 0. Classes 6, 7 and 8 
are new and not previously discussed in the literature. 

3 Renormalization group running 

The neutrino mixing parameters are measured at low energies, while the effective mass matrix 
of Eq. (1) can arise at a much high energy scale. The evolution of the light neutrino masses 
from the lightest RH neutrino mass scale Mi to the electroweak scale Mz is described by 
the one-loop RGE [13], 

IbTT^^ = kM + + M{YiYiY] , (4) 
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where t = ln(A/Mi), A is the renormalization scale and Yi = diag(|/e, I/t) with being 
the eigenvalues of the charged lepton Yukawa coupling matrix. In the SM, 
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and in the Minimal Supersymmetric Standard Model (MSSM), 
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( 6 ) 


where gi^ g 2 are the gauge couplings, yt is the top quark Yukawa coupling, and A is the Higgs 
self-coupling. 

The solution to Eq. (4) can be written in the form [14], 


M{Mz) = lo 
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where 


and 


Jo = exp 


r\n{Mi/Mz) 




K,(t)dt 


Ii = exp 




^ ' yfm 


( 8 ) 

( 9 ) 


for I = e, y, r. 

From Eq. (7), we see that RGE running can be described by multiplying the elements of 
the mass matrix by three factors, i.e.. 


May — M'^pI^Ialy , ( 10 ) 

where a, (d = e,y,T, M and M' are the effective mass matrix at the electroweak scale and 
the lightest RH neutrino mass scale, respectively. Taking the inverse of Eq. (7), we get 

(M-%f = (M'-'UjM . ( 11 ) 

Since Cay = det it is easy to verify that both the standard and nonstandard 

constraints are stable under one-loop renormalization group running from Mi to Mz- 
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Note that above the lightest RH neutrino mass, one must consider seesaw threshold 
effects, which can lead to large corrections to Qyi and the Dirac and Majorana phases [13]. 
Due to these effects, some textures that are excluded by data may be allowed above the seesaw 
threshold [15]. Estimating the effect of RGE running above Mi on low-energy parameters is 
beyond the scope of this paper because additional assumptions about the Yukawa coupling 
matrix and the RH Majorana mass matrix need to be invoked [16]. 

It has been shown that the two-loop RGE can change the rank of the light neutrino mass 
matrix [17], which implies that the constraint det M = 0 is not stable at two loop. In the 
SM, two-loop effects can be described by adding an additional term to Eq. (4) [18], so that 

16ir"^ = kM + vlW'')M + M(Y,Y*f] + M(Y,Y,')M(Y,Y<f . (12) 

The solution to the two-loop RGE is [19] 

A/(Mz) = /„ . (13) 

_M^rVrIer _ 

where 

2 j'ln{Mi/Mz) 

imn = exp -■/o yiit)yUt)dt , (i4) 

for m,n = e, /i, r. We see that the nonstandard constraints are not stable under two-loop 
running. However, as shown in Ref. [18], the two-loop correction to the lightest mass is of 
order 10“^^ eV in the SM and 10“^°(tan/5/10)"‘ eV in the MSSM. Two-loop effects can also 
generate a tiny value for ^13 at the level of 10“^^ — 10“^^ in the SM [19]. Gonsidering current 
uncertainties in the oscillation parameters, two-loop corrections do not affect our numerical 
results. 

4 Phenomenology of the nonstandard constraints 

The phenomenology of the standard constraints have been studied extensively in the lit¬ 
erature; for recent analyses, see Ref. [ 6 ] for one element zero. Ref. [7] for one off-diagonal 
element/cofactor zero, and Ref. [ 8 ] for one diagonal element/cofactor zero. The results for 
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two element/cofactor zeros can be found in Ref. [9]; other recent analyses can be found in 
Ref. [10] for two element zeros, Ref. [11] for two cofactor zeros, and Ref. [12] for one ele¬ 
ment zero and one cofactor zero. Here we focus on the phenomenology of the nonstandard 
constraints. 

4.1 Classes 1, 2, 3, 9: det M = 0 

Structures with the constraint det M = 0 are easiest to analyze. Since det M = 0 is equivalent 
to mi = 0 (m 3 = 0) for the NH (IH), the oscillation parameters are not affected, and these 
models are the same as models that have already been studied, with the additional constraint 
that the lightest mass is zero. For Class 9, det M = 0 is the only constraint and therefore any 
value of S is allowed (see Sec. 5 for a discussion of the possible constraints from neutrinoless 
double beta decay). For models with constraints besides det M = 0, we only need to check 
the parameter space allowed by the additional constraints and see if mi = 0 (m 3 = 0 ) is 
allowed for the NH (IH). 

Class 1 constraints have the same number of parameters as the most economical seesaw 
model [20]. From Ref. [9], we know that if two diagonal elements are both equal to zero, 
only = M-rr = 0 for the IH is allowed, but in this case the 2(7 allowed region does not 
include m 3 = 0. Hence, Class 1 is not allowed by the experimental data at 2a. 

For Class 2, from Ref. [8] we know that only = 0 or M^- = 0 are allowed for the 
IH at 2a when the lightest mass is zero. Similarly, for Class 3, from Ref. [7] we know that 
only Mefx = 0 or M^t = 0 are allowed for the IH at 2a when the lightest mass is zero. In 
all other cases with exactly one texture zero in the light neutrino mass matrix, a massless 
neutrino is excluded at 2a. Hence only four cases in Class 2 and 3 are allowed for the IH at 
2 ( 7 , and none for the NH. The allowed cases in Class 2 are: (i) = 0, and det M = 0; (ii) 

Mrr = 0, and det M = 0. The allowed cases in Class 3 are: (i) = 0, and det M = 0; (ii) 

Mer = 0, and detM = 0. We note that the Class 2 and 3 constraints have been noted in 
Ref. [4] as Eqs. (48) and (49). Our results for the allowed cases are consistent with Ref. [4]. 

Reference [4] also lists the conditions det M = 0 and — {Ma/sY = 0 (as Eq. 50) 

and Maa = = 0 and {M~^)aa = 0 (as Eq. 51) as new constraints. The structures that 

lead to these conditions have six and three cofactor zeros, respectively. We do not count 



them as nonstandard constraints because a light neutrino mass matrix with more than two 
cofactor zeros has already been shown to be excluded by experimental data. 


4.2 Class 4 


Our Class 4 constraints are equivalent to M^a = 0, = 1 and ^ 0 

obtained in Ref. [4] as Eq. (51). We demonstrate this equivalency below. 

Since M^a = 0, 

det M = . (15) 

Plugging MpisMa-f = into the above equation, we get 

deiM = (16) 

Also, since Maa = 0, we have , and therefore 

= 1. (17) 

Together with M^a = 0 and ^ 0 (which is implicitly included in our constraints), we 
see Class 4 constraints are consistent with the conditions (52) in Ref. [4]. The form of the 
constraint in Table 1 has the merit of being quadratic instead of cubic (as in Ref. [4]), which 
allows for easier numerical analysis. 

The constraints M^a = 0 and can be written as 

Xi + < 7 X 2 + pX^ = 0 , (18) 


and 


Ai + cA'A 2 + A^ + < 7 Ai2 + pAi3 + crpA23 — 0 , (19) 

where a = (m2/mi)e"*‘^^ p = (m3/mi)e“*'^h Xi = Ai = Up^U*iU*i and Aij = 

2U;,U*^^iU*,U;^+U*^U;,)-^^^^^ for * = l, 2 , 3 . Since Eq. (19) is quadratic, 

it can be solved with Eq. (18). We get 


423 ^ 1 ^ 2^3 - A13XIX3 - 2A2XyXi + A12X2XI ± X2 
2X3(43X1 + X 3 ( 42 A '3 — ^23^2)] 

423 ^ 1^3 + 413X2X3 — 243X1X2 — 4 i 2 X| =P \/A 
2(43X1 + X 3 ( 42 X 3 — 423X2)] 


( 20 ) 

( 21 ) 
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nil (eV) m3 (eV) 

Figure 1: The 2a allowed regions in the plane of the lightest mass and Dirac phase for the 
constraints, = 0 and = 2M^T-Me,-. The left (right) panel is for the normal 

(inverted) hierarchy. 

where A = [2A^XiX2+X^{Ai2X^ — A 2 ^Xi — Ai^X 2 )Y —^[A^Xl+X^{AiX^ — Ai^Xi)\[A^X'^ + 
X 3 {A 2 Xs — A 23 X 2 )]. Taking the absolute values of a and p, we hnd m 2 /mi and mz/rrii. 
Then mi can be written as 



where the plus (minus) sign in Eq. (22) is for the NH (IH). 

Our numerical analysis proceeds as follows. Since there are 5 independent parameters 
in the light mass matrix in this case, we take them to be 6 ^ 12 , ^ 23 , 6 ^ 13 , 5m^ = and 

Xm? = m\ — {m\ + m 2 )/ 2 . First, we choose a set of the hve oscillation parameters within 
the 2a range from a recent global £t [21]. Then we solve Eq. (22) by scanning over 5 from 0 
to 360°. We sample 10^ sets of oscillation parameters allowed at 2a] by keeping the values 
of mi and 5 that are allowed, we obtain the allowed regions in mi — 5 (m 3 — 5) plane for the 
NH (IH). 

The allowed cases in Class 4 are as follows: 

1. = 0 and IVlrrM^^ = 2M^,-Mer. This case is allowed for both hierarchies. The 
allowed regions are shown in Fig. 1. 

2. Mrr = 0 and = 2M^,-Me^. This case is allowed for both hierarchies. The 

allowed regions are shown in Fig. 2. They are similar to Fig. 1 with 5 —)■ 5 + 180° for 
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nil (eV) m3 (eV) 

Figure 2: Same as Fig. 1, except for M^r = 0 and M^^Mer = 2M^T-Me^. 

both hierarchies. The similarity is explained by the approximate fi — t symmetry of 
the experimental data. 

Note that these allowed regions are consistent with the results obtained in Ref. [4], in 
which 3a ranges of each parameter are used as input. 

4.3 Class 5 

Similar to Class 4, Class 5 constraints are consistent with the conditions {M~^)aa = 0, 
Mi 3 is{M~^)i 3 i 3 = 1 and 7 ^ 0 in Eq. (53) of Ref. [4]. We hnd that it is easier to use the 
quadratic form of Class 5 constraints for the numerical analysis. 

Similar to Eqs. (18) and (19), the two constraints Caa = 0 and = 2Cq,^Cq,..), can 

be written as 

Xi + X2/a* + X3/p* = 0, (23) 

and 

Ai + A 2 /{a*)'^ + 43 /(p*)^ + Ayija* + Ais/ p* + R 23 /(<^*P*) = 0 ^ (24) 

where p, a, Xj, A^ and Aij are dehned after Eq. (19). The analysis of this case follows that 
of Class 4 with the replacements, a —)■ 1/a* and p —)■ 1/p*. The allowed cases in Class 5 are 
as follows: 

1. Cee = 0 and C^^Cer = 2Ce^C^r- This case is allowed for the IH only. The allowed 
regions are shown in the left panel of Fig. 3. 
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Figure 3: The 2a allowed regions in the { 7713 , 6 ) plane for the inverted hierarchy. The 
left (right) panel shows the results for the constraints, C^e = 0 and 



m\ (eV) m3 (eV) 


Figure 4: Same as Fig. 1, except for = 0 and CrrCe^ = 2C^rCeT- 

2. Cee = 0 and CrrCen = 2CerC'^T- This case is allowed for the IH only. The allowed 
regions are shown in the right panel of Fig. 3, which is similar to the left panel with 
(5 —)■ (5 + 180° because of the approximate /i — r symmetry. 

3- = 0 and CrrCe^ = 2C^rCeT- This case is allowed for both hierarchies. The allowed 

regions are shown in Fig. 4. 

4. Crr = 0 and C^^Cer = 2Cf^T-Ceii- This case is allowed for both hierarchies. The allowed 
regions are similar to Fig. 4 with 5 —)■ 5 + 180° for both hierarchies. The similarity can 
be explained by the approximate /i — r symmetry of the experimental data. 

Note that these allowed regions are consistent with the results obtained in Ref. [4]. Also, 
the Class 4 constraints are dual to the Class 5 constraints because they have the same 
functional form with the roles of elements and cofactors reversed [22]. From Ref. [22], we 
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know that the allowed regions of Class 4 and 5 should be similar for opposite mass hierarchies 
when the lightest mass is larger than about 20 meV. This can be seen from the left panel of 
Fig. 1 and the right panel of Fig. 4, and from the right panel of Fig. 1 and the left panel of 
Fig. 4. 

4.4 Class 6 

Since detM ^ 0, the constraint MaaCaa = detM is equivalent to Maa{M~^)aa = 1, which 
can be written as 


(miC:? + = 1 • (25) 

Dehning Uaj = \Uaj\e^^^ and 0' = (f)j + 2/3j — 2[3i , we absorb the phases of the f/’s into (p/s 
and write Eq. (25) as 

(m,|£;„,p + m 2 |C/„ 2 pe""*« + m 3 |C/„ 3 pe""*^(m;''|£;„, p + + m. 3 'IC/^sl V*'“) = 1 

(26) 

We now expand the above equation into real and imaginary parts, and after some simplih- 
cation we get 


\Ual\^ + \Ua2\^ + \Ua3\* + \Ua2\‘^\Ua3\^ ( — + —^ (COS 02 COS 03 + siu 02 siu 03) 

' “l3 7712/ 


+ \U.l\^\U^2\^ + COS0'2 + |f/«ir|f/„3r COS0' = 1 , (27) 

\m2 rriij yms mi J 


for the real part and 


irr |2|rr \2 f 7712 \ . , . |2|_ |2/^"^1 "^3 \ . ,/ 

\Ual\ \Ua2\ -Sm02 + \Ual\ \Ua3\ -SIU 03 

\7772 7771 / \7773 7771 / 

+ \Ua2\‘^\Ua3\‘^ f— “ —^ (cOS 03 slu 02 - sin 03 COS 02 ) = 0, (28) 

07772 7773 / 


for the imaginary part. The imaginary part can be put into the form 


D = A cos 02 + B sin 02 , 


(29) 
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where 


^ Sin 0 ' 3 , 

\m2 msj 

U-\TT |2|rr |2/^"^3 m2\ , I | 2 |.. , 2(^1 "^2 A 

B — \Ua2\ lUaSl I-COS03 + \Ual\ \Ua2\ I-, 

\m2 m3 J \m2 mi J 

n ITT |2|rr 12 . ,/ 

D = -\Ual\ \Ua3\ - 811103 . 

\m3 mij 

Equation (29) has the solution (as long as 


sin 02 
cos 02 


BD ± + B^-D^ 

^2 + ^2 

adt 

A2 + fi2 


(30) 

(31) 


For a given 03 we solve for 02 and plug that into Eq. (27). Then for a hxed set of S, 
mi {m 3 ) for the NH (IH), and a given set of oscillation parameters, we scan 03 to see where 
(if) Eq. (27) is satished. We also ignore those values of 03 that have no real solution for 
00 i.e., A'^ + B‘^ < D^. If a solution is found, that point in mi — 6 space is allowed for the 
corresponding set of oscillation parameters. With the value of 03 that gives the solution, 02 
and 03 are also obtained. 

The allowed cases in Class 6 are as follows: 


1. MeeCee = det M. This case is allowed for both hierarchies. There is no constraint on 
5 in this case. For the best-£t oscillation parameters, the allowed range of the lightest 
mass is mi > 0.0024 eV for the NH and m 3 > 0.00050 eV for the IH. The allowed 
range at 2(j is mi > 0.0020 eV for the NH and m 3 > 0.00037 eV for the IH. 

2. = detM. This case is allowed for both hierarchies, and the allowed regions 
are shown in Fig. 5. 

3. MrrCrr = det M. This case is allowed for both hierarchies, and the allowed regions 
are shown in Fig. 6 . The allowed regions are similar to Fig. 5 with 5 —)■ 5 + 180° for 
both hierarchies because of the approximate /i — r symmetry of the data. 
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nil (eV) m3 (eV) 

Figure 5: The allowed regions for the constraint = detM. The left (right) panel 

is for the normal (inverted) hierarchy. The dark shaded regions correspond to the best-fit 
oscillation parameters, while the light shaded regions are allowed at 2a. 



mi (eV) m3 (eV) 

Figure 6: Same as Fig. 5, except for MrrCrr = det M. 


15 






4.5 Class 7 


With X = and y = the constraint M'^^Caa = can be written as 

ax^ + bx'^ + ex + d = 0 , (32) 

where 

™2rr* 2 

“ " ^ - 4(mil7^3C',s + , (33) 

l> = {2m.(mX^3 + ^mUl,v)U-jK, 

rriim^ 

-AmiimiUpzU^z + m‘iUyiU^iy)Ul^{Up2U*^ + Ul^U*2) 

+m,y [m,Ul^U*J + “ 2 m,U^ 2 U, 2 U;,U;,) 

-AiTmUysU.s + msyUf,,U,,)U:,iU;,U;, + U;,U;,)] } , (34) 


c = 


rriimz 


{ml{m,Ul^ + m,Ul,y)Wj 


+2mi[/*i^ [m3|/(mi|[/„2|^ + mi\Uaz\^ + msyUl^U*^"^) 
-2mi{miU/33U^3 + m3yUyiU^i)U^^U*^] 

- Am.m^yU:, [m,Up2U,2K,{U*p^U;, + U;,U;^) 

HmiUf^sU.s + m3yUf,,u,,)u:3iu;3u;, + u^.u;^)] 

+mly^U:, [2m,\U^2\^U:, - 4{m,Up3U,3 + m3yUisiU^,)U:,U;,U;, 

HmiUl, + m3Ul,y)iu:,f - + u;,u;,)]} , 


(35) 


d = 


JL 

m 2 


[^lUl2U*J + 2m,U*J{m3Ul,yU*J - 2m,Uy2U,2U;,U^ 


7I J 


-im,m.,UfX-,2yU:Xi>(UhU;i + UyX;^) + mWU'JXlX'J - iU„X.flUkU;^)\ 


(36) 


For a hxed set of 5, mi (m 3 ) for the NH (IH), and a given set of oscillation parameters, 
we solve the cubic equation (Eq. 32) to hnd x as a function of y, and then scan over 03 from 
0 to 27r to hnd where (if) |x| = 1. If |x| = 1 is found, there is a solution for 02 and 03. 
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By keeping the values of mi (m 3 ) and S that have a solution, we obtain the allowed regions 
in the mi — 6 (m 3 — 5) space for the NH (IH). All the cases in Class 7 are allowed for the 
both hierarchies, and the allowed regions are shown in Fig. 7. The 2a allowed regions in 
the second row of Fig. 7 are similar to the third row of Fig. 7 with <5 —)■ 5 + 180° for both 
hierarchies because of the approximate /i — r symmetry of the data. 


4.6 Class 8 

Since det M 7 ^ 0, the constraint det M = —M^i^Caa is equivalent to 

Similar to Class 7, dehning x = and y = we can write the constraint as 

a'x^ + h'x"^ + c'x + d' = 0 , 


where 


a =- 

minis 


(37) 


(38) 


(39) 


b' = [vu;, (|C/„2| ‘ + 1) + 2u:, (!/mr‘£7"i + "> 3 " imu:xh + ’^3VU:,Uk)] , 

(40) 


+ (™i£7olC7p + rn3yU‘„,UUf + ™3‘C^os) . 


(41) 


<i' = 


m 2 


(42) 


The solutions for 02 and 03 then proceed as in Class 7. We hnd all the cases in Class 8 
are allowed for the both hierarchies and the allowed regions are shown in Figs. 8 and 9. Due 
to the approximate y — r symmetry in the experimental data, the 2a allowed regions in each 
row of Fig. 8 are similar to those in the corresponding row of Fig. 9 with <5 —)■ <5 + 180° for 
both hierarchies. 
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Figure 7: The first, second and third rows show the allowed regions for the constraints, 
M^^Cee = 4:Mef,MerC^r, = AM^^M^^Cer and M^^Crr = AMerM^^Cef,, respectively. 

The dark shaded regions correspond to the best-£t oscillation parameters, while the light 
shaded regions are allowed at 2(j. The left (right) panels are for the normal (inverted) 
hierarchy. 
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Figure 8: Same as Fig. 7, except that the hrst, second and third rows show the results for the 
constraints, M^^detM = —M^^Cee, MeedetM = —and MrrdetM = 
respectively. 


19 
























Figure 9: Same as Fig. 7, except that the hrst, second and third rows show the results for 
the constraints, M^-j. det M = —M^^Cee, M^e det M = —M^^Crr and det M = 
respectively. 
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5 Discussion and conclusion 


We carried out a systematic study of zeros in the Dirac and RH Majorana neutrino mass 
matrices in the context of the Type I seesaw mechanism. We derived complex constraints 
on the light neutrino mass matrix for various textures in both the Dirac and RH Majorana 
neutrino mass matrices. There are 9 nonstandard constraints besides the standard ones 
that can be expressed in the form of one or two element/cofactor zeros alone. We showed 
that both the standard and nonstandard constraints are stable under one-loop RGE running 
from the lightest RH neutrino mass scale Mi to the electroweak scale Mz- In addition, we 
studied the phenomenological implications for the nonstandard constraints, and found that 
some cases for the normal or inverted hierarchy are excluded, and for the rest we obtain 
the allowed regions for the lightest mass and Dirac CP phase, which will be probed in the 
next generation of neutrino experiments. We found 12 new models (Classes 6, 7 and 8) 
not previously discussed in the literature which are allowed at 2a for both the normal and 
inverted hierarchies. 

Once the lightest neutrino mass and Dirac CP phase are determined, the rate for neutri¬ 
noless double beta decay (0z//9/9), which depends on the magnitude of the z/g — r'e element of 
the neutrino mass matrix, 

|Mee| = \mic\2c\^ + , (43) 

is hxed. 

In Table 2, we list the minimum and maximum values of |Mee| at the 2a level for Class 2, 3, 
4 and 5 constraints. In Table 3, we list the minimum values for Class 6, 7, 8 and 9 constraints 
with the best-£t oscillation parameters and the 2a lower bounds. The maximum values for 
Class 6, 7, 8 and 9 constraints are all above 1000 meV. The most stringent experimental 
upper limit on the effective mass \Mee\ is 120 — 250 meV at 90% C.L. [23]. Some cases in 
Class 3, 4 and 5 could be ruled out once the sensitivities to |Mee| reach about 50 meV for 
the currently running 0z//5/9 experiments [24]. Nevertheless, the minimum values of |Mee| for 
Class 6, 7, 8 and 9 constraints are all below 50 meV and can be completely probed only if 
the sensitivity of future 0z//9/9 experiments is significantly improved. Note that the minimum 
value of |Mee| is about 15 meV for the IH and 0 for the NH. The mass hierarchy will also be 
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Class 

Constraints 

Hierarchy 

Minimum 

Maximum 

2 

= 0 and det M = 0 

IH 

15 

18 

2 

Mtt = 0 and det M = 0 

IH 

15 

17 

3 

Me^ = 0 and det M = 0 

IH 

46 

48 

3 

Mer = 0 and det M = 0 

IH 

46 

48 

4 

= 0 and = 2M^rMer 

NH 

41 

278 

4 

= 0 and M^rMe^ = 2M^rMeT 

IH 

16 

214 

4 

Mrr = 0 and 

NH 

65 

277 

4 

Mrr = 0 and 

IH 

15 

223 

5 

Cee = 0 and C^^Cer = 2CeiiCftr 

IH 

19 

23 

5 

Cee = 0 and CrrCeii = 2(CerC^r 

IH 

19 

24 

5 

C^^ = 0 and CrrCe^ = 2C^rCeT 

NH 

64 

278 

5 

Cfifi = 0 and CrrCe^ = 2C^rCeT 

IH 

65 

222 

5 

Ctt = 0 and C^^Cer = 2(C^T-Ce^ 

NH 

39 

268 

5 

Ctt = 0 and C^^Cer = 2(C^T-Ce^ 

IH 

61 

219 


Table 2: The minimum and maximum values of |Mee| (in meV) for Class 2, 3, 4 and 5 
constraints at 2a. 

determined if |Mee| can be pushed to 15 meV in the future. 

Recent global hts [21, 25, 26] indicate a slight preference for a nonzero Dirac CP phase, 
with a preferred value S ~ 270°. Also, combined T2K and reactor analyses show a weak 
hint for the normal hierarchy [27]. Most viable cases of the nonstandard constraints allow 
6 = 270°, except for two cases in Class 5: Cee = 0, C^^Cer = 2Ce^C^r and Cge = 0, 
CrrCe^ = 2CeT-C^T-. These two cases may be ruled out if future measurements hnd 5 to be 
close to 270°. If the hint for the normal hierarchy is conhrmed. Classes 2 and 3, and the IH 
cases in the other classes will be ruled out. 

Acknowledgements. This research was supported by the DOE under grant No. DE- 
SC0010504 and by the Kavli Institute for Theoretical Physics under NSF grant No. PHYll- 
25915. 
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Class 


Constraints 

Best-fit 

2a lower bound 

NH 

IH 

NH 

IH 

MeeCee = det M 

3.2 

17.9 

2.6 

15.1 

= det M 

0.2 

17.9 

0.0 

14.9 

Mt-tCtt = det M 

0.8 

17.9 

0.0 

15.0 

MIC,, = 4M,^M,rC^r 

0.0 

41.4 

0.0 

39.6 

= 4M^rM,^C,r 

0.9 

17.9 

0.6 

15.0 

M^Crr = 4M,rM^^C,^ 

0.8 

17.9 

0.5 

15.0 

det M = 

1.6 

18.0 

1.3 

14.9 

Mrr det M = -MlC„ 

1.7 

17.9 

1.4 

14.9 

M„detM = 

0.0 

21.7 

0.0 

18.1 

MrrdeiM = -Ml^C^^ 

0.5 

17.9 

0.2 

15.0 

M„ det M = -M^Crr 

0.0 

20.3 

0.0 

18.6 

= -Ml^Crr 

0.5 

17.9 

0.2 

15.0 

det M = 0 

1.4 

17.9 

1.2 

14.8 


Table 3: The minimum values of |Mee| (in meV) for Class 6, 7, 8 and 9 constraints for the 
best-fit oscillation parameters, and the 2a lower bounds. 
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A Textures for the nonstandard constraints 


We list all possible structures that lead to the nonstandard constraints for N = + > 7. 

Also we do not consider > 7 because it has one row of zeros in M^), which leads to three 
zeros in M, which is excluded by experimental data. Following Ref. [9], we dehne Class X 
constraints on M or C = det(M)M“^ as one zero on diagonal and one off-diagonal zero 
sharing column and row, Class Y constraints as one zero on diagonal and one off-diagonal 
zero not sharing column and row, and Class Z constraints as two zeros on diagonal. 


A.l Nfi — 4 


For Nr = 4, there is only one structure for Mr, 


X 


Mr = 


0 

0 


0 0 
0 X , 
X 0 


(44) 


where the symbol x denotes a nonzero matrix element. All other structures of Mr can be 
obtained by a permutation of the columns and rows of the above structure. If Mr has more 
than 5 zeros, M or C would have at least three zeros, so here we only consider Nr < 5. 


1. Nr = 5. Of the 126 different structures of Mr in this case, 120 lead to at least three 
zeros in M or C, and only 6 structures lead to Class 1 nonstandard constraints. The 
6 structures are: 



0 

X 

0 


0 

0 

X 

Mr = 

0 

0 

X 


0 

X 

0 


0 

X 

X 


0 

X 

X 


which lead to the constraints, M^e = = 0 and det M = 0; 



0 

X 

0 


0 

0 

X 

Mr = 

0 

X 

X 


0 

X 

X 


0 

0 

X 


0 

X 

0 


(45) 


( 46 ) 


24 



which lead to the constraints, Mge = = 0 and det M = 0; 



1 

X 

X 

O 


1 - 

X 

X 

O 

Mo = 

0 x0 


X 

O 

O 


0 0 X 


-1 

o 

X 

O 


which lead to the constraints, = M^-t- = 0 and det M = 0. 


(47) 


2 . No = 4. Of the 126 different structures of Mo in this case, 39 lead to at least three 
zeros in M or C, 3 lead to a block diagonal matrix for M, 36 lead to Class X constraints 
on M, 6 lead to Class Z constraints on M, 6 lead to Class Z constraints on C, and 18 
lead to one zero in M and one zero in C. In addition, there are 6 structures of Mo 
that lead to Class 2 nonstandard constraints: 



0 

X 

0 


0 

0 

X 

Mo = 

0 

X 

X 


0 

X 

X 


0 

X 

X 


0 

X 

X 

which lead to the constraints, Mee 

= 

0 and det M 

= 

0 ; 



0 

X 

X 


0 

X 

X 

Mo = 

0 

X 

0 


0 

0 

X 


0 

X 

X 


0 

X 

X 

which lead to the constraints, M^^ 

= 

0 and det M = 

0 ; 



0 

X 

X 


0 

X 

X 

Mo = 

0 

X 

X 


0 

X 

X 


0 

X 

0 


0 

0 

X 


which lead to the constraints, Mrr = 0 and det M = 0. 


(48) 


(49) 


(50) 


The remaining 12 structures lead to Class 4 nonstandard constraints: 



0 

X 

0 


0 

0 

X 

to 

0 

X 

X 


0 

X 

X 


X 

0 

X 


X 

X 

0 


(51) 
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which lead to the constraints, 


Mee = 0 and 




0x0 

Md = X 0 X , 


0 X 
X 0 , 


0 X X 


which lead to the constraints, 


Mee = 0 and MrrMe^ 


2MerM^r; 


0 X X 

Md = 0 X 0 , 


which lead to the constraints. 


X 0 X 


= 0 and 


X 0 X 


Mo = 0 X 0 , 

0 X X 


X X 


0 X , 

X 0 

= 2Me^Me,; 

X 0 
0 X , 


X X 


which lead to the constraints. 


M^^ = 0 and 


= 2M.Mer] 


X 0 X 

Mo = 0 X X , 

0x0 


X 0 


X X 


0 X 


which lead to the constraints. 


Mrr = 0 and M^^Mer 


= 2Me^M^,; 


0 X X 

Mo = X 0 X , 

0x0 


X 0 , 

0 X 


which lead to the constraints. 


Mrr = 0 and MgeM^T- 


2Me^Mer. 


3. Nd = 3 . Of the 84 different structnres of Mo in this case, 5 lead to at least three 
zeros in M or O, 12 lead to Class X constraints on M, 24 lead to one diagonal zero 
in M, 6 lead to one off-diagonal zero in M and 24 lead to one diagonal zero in C. In 
addition, there are 6 structnres that lead to Class 6 nonstandard constraints: 


X X 0 

Mo = 0 X X , 


X 0 X 

0 X X , 


0 X X 


0 X X 
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which lead to the constraint, M^eCee = detM; 



0 

X 

X 


0 

X 

X 

Md = 

X 

X 

0 


X 

0 

X 


0 

X 

X 


0 

X 

X 


which lead to the constraint, = det M; 



0 

X 

X 


0 

X 

X 

Mo = 

0 

X 

X 


0 

X 

X 


X 

X 

0 


X 

0 

X 


which lead to the constraint, M^T-Crr = det M. 

Also, there are 6 structures lead to Class 7 nonstandard constraints: 



1 

O 

X 

X 

_ 1 


1- 

X 

X 

O 

Md = 

o 

X 

X 


X 

o 

X 


1 

X 

O 

X 

_1 


-1 

O 

X 

X 

_1 


which lead to the constraint, 



X 

X 

0 


X 

0 

X 

Md = 

0 

X 

X 


0 

X 

X 


X 

0 

X 


X 

X 

0 


which lead to the constraint, = iMefiM^rCer] 



i 

o 

X 

X 


1- 

X 

o 

X 

Md = 

X 

o 

X 


o 

X 

X 


1 

o 

X 

X 

1 _ 


1 

o 

X 

X 

1 _ 


which lead to the constraint, M^^Crr = dMerM^^Ce^. 

The remaining one structure leads to the Class 9 constraint, det M = 0: 


( 58 ) 


(59) 


(60) 


(61) 


(62) 


0 

X 

X 

0 

X 

X 

0 

X 

X 


( 63 ) 
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A.2 Nr = 3 


For Nfi = 3, there are four structures of that are not equivalent under permutation. 
They are 


Mri = 


Mrs — 


X 

0 

0 


X 

0 

X 

0 

X 

X 

, Mr2 — 

0 

0 

X 

0 

X 

0 


X 

X 

0 

0 

X 

X 


X 

0 

0 

X 

0 

X 

, Mr4 = 

0 

X 

0 

X 

X 

0 


0 

0 

X 


(64) 


All other structures of Mr can be obtained by a permutation of the columns and rows of 
the above structure. If Mr has more than 5 zeros, M or C would have at least three zeros 
for all three structures of Mr, so here we only consider Nr < 5. 


A.2.1 Mri 

1. Nr = 5. Of the 126 different structures of Mr in this case, 84 lead to at least 


three zeros in M or C, 6 lead to a block diagonal matrix for M, 12 lead to Class X 
constraints on M, and 12 lead to one zero in M and one zero in C. In addition, there 
are 6 structures that lead to Class 2 nonstandard constraints: 



0 

X 

0 


0 

X 

0 


Mr = 

0 

X 

X 


0 

0 

X 

, (65) 


0 

0 

X 


0 

X 

X 


which lead to the constraints, Mge 

= 

0 and det M 

= 

0; 




0 

X 

X 


0 

0 

X 


Mr = 

0 

X 

0 


0 

X 

0 

, (66) 


0 

0 

X 


0 

X 

X 


which lead to the constraints, M^^ 

= 

0 and det M = 

0; 




0 

0 

X 


0 

X 

X 


Mr = 

0 

X 

X 


0 

0 

X 

(67) 


0 

X 

0 


0 

X 

0 
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which lead to the constraints, M^-j. = 0 and det M = 0. 

The remaining 6 strnctnres lead to Class 3 nonstandard constraints: 



X 

0 

0 


0 

0 

X 

Mo = 

0 

0 

X 


X 

0 

0 


X 

0 

X 


X 

0 

X 


which lead to the constraints, = 0 and det M = 0; 



X 

0 

0 


0 

0 

X 

Mo = 

X 

0 

X 


X 

0 

X 


0 

0 

X 


X 

0 

0 


which lead to the constraints, M^r = 0 and det M = 0; 



X 

0 

X 


X 

0 

X 

Mo = 

X 

0 

0 


0 

0 

X 


0 

0 

X 


X 

0 

0 


which lead to the constraints, = 0 and det M = 0. 


( 68 ) 


(69) 


(70) 


2. No = 4. Of the 126 different strnctnres of Mo in this case, 33 lead to at least three 
zeros in M or C, 3 lead to a block diagonal matrix for M, 18 lead to Class X constraints 
on M, and 6 lead to one zero in M and one zero in C, 15 lead to one diagonal zero in 
M, 12 lead to one off-diagonal zero in M, 9 lead to one diagonal zero in C, and 6 lead 
to one off-diagonal zero in C. In addition, there are 3 strnctnres that lead to Class 2 
nonstandard constraints: 


Md = 


0 

0 

0 


X 

X 

X 


0 

X , 
X 


(71) 


which leads to the constraints, Mee = 0 and det M = 0; 


Mo = 


0 X X 
0x0 
0 X X 


(72) 
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which leads to the constraints, = 0 and det M = 0; 

0 X X 

Mo = 0 X X , (73) 

0x0 

which leads to the constraints, M^-t = 0 and detM = 0. 

Also, there are 6 structures that lead to Class 6 nonstandard constraints: 

X X 0 X X 0 

Md= 0 X X , OOx, (74) 

0 0 X 0 X X 

which lead to the constraint, M^eCee = detM; 

0 X X OOx 

Mo = X X 0 , X X 0 , (75) 

OOx 0 X X 

which lead to the constraint, = det M; 

0 X X OOx 

Mo = 0 0 X , 0 X X , (76) 

X X 0 X X 0 

which lead to the constraint, MrrCrr = det M. 

There are 6 structures lead to Class 8 nonstandard constraints: 

X X 0 

Mo = 0 0 X , (77) 

X 0 X 

which leads to the constraint, M^^detM = —M^^Cee] 

X X 0 

Mo = X 0 X , (78) 

OOx 
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which leads to the constraint, Mn- det M = —M^^Cee] 

0 0 X 

Md = X X 0 , 

X 0 X 

which leads to the constraint, MeedetM = 

X 0 X 

Md = X X 0 , 

0 0 X 

which leads to the constraint, Mrr det M = — 

X 0 X 

Mb = 0 0 X , 

X X 0 

which leads to the constraint, M^^detM = 

0 0 X 

Mb = X 0 X , 

X X 0 

which leads to the constraint, MeedetM = —M^^Crr- 


The remaining 9 structures lead to the Class 9 constraint, det M = 0: 


X 0 
Mb = X 0 
X 0 

X 0 
X 0 
0 0 

0 X 
0 X 
0 0 



X 0 X 

X 0 0 , 

X 0 X 

X 0 X 

0 0 X , 

X 0 X 

0 X X 

0 0 X , 

0 X X 


X 

0 

0 

X 

0 

X 

X 

0 

X 

0 

0 

X 

X 

0 

X 

X 

0 

X 

0 

0 

X 

0 

X 

X 

0 

X 

X 


(79) 


(80) 


(81) 


(82) 


(83) 


(84) 


(85) 
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A .2.2 

1. Nd = 5. Of the 126 different structures of in this case, 90 lead to at least three 
zeros in M or C, 12 lead to Class X constraints on M, 6 lead to Class Z constraints on 
C, and 6 lead to one zero in M and one zero in C. In addition, there are 6 structures 
that lead to Class 2 nonstandard constraints: 



0 

0 

X 


0 

0 

X 

Md = 

0 

X 

0 


0 

X 

X 


0 

X 

X 


0 

X 

0 


which lead to the constraints, Mee = 0 and det M = 0; 



0 

X 

0 


0 

X 

X 

Md = 

0 

0 

X 


0 

0 

X 


0 

X 

X 


0 

X 

0 


which lead to the constraints, = 0 and det M = 0; 



0 

X 

0 


0 

X 

X 

Md = 

0 

X 

X 


0 

X 

0 


0 

0 

X 


0 

0 

X 


which lead to the constraints, = 0 and det M = 0. 


( 86 ) 


(87) 


( 88 ) 


The remaining 6 structures lead to Class 5 nonstandard constraints: 


Md = 


0 X X 
0x0 
X 0 0 


which leads to the constraints, C^e = 0 and = 2Ce^C^r, 


Md = 


0 X X 
X 0 0 
0x0 


(89) 


(90) 
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which leads to the constraints, Cee = 0 and CrrCefj, = 2CerC'/iT; 

0x0 

Md = 0 X X , (91) 

X 0 0 

which leads to the constraints, = 0 and CeeC^r = ‘^Ce^Cer] 

X 0 0 

Md = 0 X X , (92) 

0x0 

which leads to the constraints, = 0 and CrrCgf^ = ^C^rCer'^ 

X 0 0 

Mb = 0 X 0 , (93) 

0 X X 

which leads to the constraints, Crr = 0 and C^^Cer = ‘^C^rCe^] 

0x0 

Mb = X 0 0 , (94) 

0 X X 

which leads to the constraints, Crr = 0 and CeeC^r = ‘^CerCe^i- 

2. Nd = 4. Of the 126 different structures of Mb in this case, 39 lead to at least three 
zeros in M or O, 6 lead to Class Z constraints in C, 18 lead to Class X constraints on 
M, and 3 lead to one zero in M and one zero in C, 12 lead to one diagonal zero in M, 
and 33 lead to one diagonal zero in C. In addition, there are 3 structures that lead to 
Class 2 nonstandard constraints: 

0 0 X 

Mb = 0 X X , (95) 

0 X X 

which leads to the constraints, Mee = 0 and det M = 0; 

0 X X 

Mb = 0 0 X , (96) 

0 X X 
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which leads to the constraints, = 0 and det M = 0; 

0 X X 

Mo = 0 X X , (97) 

0 0 X 

which leads to the constraints, Mrr = 0 and detM = 0. 

Also, there are 9 structures leads to Class 6 nonstandard constraints: 

xOx xOx xOO 

Mo =0xx, 0x0, Oxx, (98) 

0x0 Oxx Oxx 

which lead to the constraint, M^eCee = det M; 

Oxx 0x0 Oxx 

Mo =x0x, xOx, xOO, (99) 

0x0 Oxx Oxx 

which lead to the constraint, = detM; 

Oxx 0x0 Oxx 

Mo =0x0, Oxx, Oxx, (100) 

xOx xOx xOO 

which lead to the constraint, MrrCrr = detM. 

The remaining 3 structures lead to the Class 9 constraint, det M = 0: 

Oxx Oxx 0x0 

Mo =0xx, 0x0, Oxx- (101) 

0x0 Oxx Oxx 

A.2.3 Mrs 

1. No = 5. Of the 126 different structures of Mo in this case, 90 lead to at least three 
zeros in M or C and 36 lead to Class Z constraints on C. 

2. No = 4. Of the 126 different structures of Mo in this case, 45 lead to at least three 
zeros in M or C, 18 lead to Class Z constraints on C, and 63 lead to one diagonal zero 
in C. 
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A.2.4 Mra 

1. Nr = 5. Of the 126 different structures of Mr in this case, 72 lead to at least three 
zeros in M or C* and 36 lead to a block diagonal matrix for M. The remaining 18 
structures lead to Class 3 nonstandard constraints: 


Md = 


which lead to the constraints, Men = 0 and det M = 0; 


Mr = 


which lead to the constraints, Mer = 0 and det M = 0; 


X 

0 

0 


X 

0 

0 


0 

X 

0 

0 

X 

0 


0 

0 

X 


0 

0 

X 

X 

X 

0 


X 

0 

X 


0 

X 

X 

0 

X 

0 


0 

0 

X 


0 

0 

X 

X 

0 

0 


X 

0 

0 


0 

X 

0 

X 

X 

0 


X 

0 

X 


0 

X 

X 


X 

0 

0 


X 

0 

0 


0 

X 

0 

X 

X 

0 


X 

0 

X 


0 

X 

X 

0 

X 

0 


0 

0 

X 


0 

0 

X 

0 

X 

0 


0 

0 

X 


0 

0 

X 

X 

X 

0 


X 

0 

X 


0 

X 

X 

X 

0 

0 


X 

0 

0 


0 

X 

0 


( 102 ) 


(103) 


(104) 


(105) 


Mo = 


X 

X 

0 


X 

0 

X 


0 

X 

X 

0 

X 

0 


0 

0 

X 


0 

0 

X 

X 

0 

0 


X 

0 

0 


0 

X 

0 

X 

X 

0 


X 

0 

X 


0 

X 

X 

X 

0 

0 


X 

0 

0 


0 

X 

0 

0 

X 

0 


0 

0 

X 


0 

0 

X 


(106) 


(107) 


which lead to the constraints, M^r = 0 and det M = 0 


2. Nr = 4. Of the 126 different structures of Mr in this case, 27 lead to at least three 
zeros in M or C, 9 lead to a block diagonal matrix for M, 54 lead to one off-diagonal 
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zero in M, and 18 lead to one off-diagonal zero in C. The remaining 18 structures lead 
to the Class 9 constraint, det M = 0: 



X 0 

X 0 , 

0 0 

X 0 

X 0 , 

X 0 

0 X 

0 X 

0 0 

0 X 

0 X , 

0 X 

X X 

X X , 

X 0 

X X 

X X , 

0 X 



X X 0 

X 0 0 , 

X X 0 

X X 0 

0x0, 

X X 0 

X 0 X 

X 0 0 , 

X 0 X 

X 0 X 

0 0 X , 

X 0 X 

0 X X 

0x0, 

0 X X 

0 X X 

0 0 X , 

0 X X 


X 

0 

0 

X 

X 

0 

X 

X 

0 

0 

X 

0 

X 

X 

0 

X 

X 

0 


X 

0 

0 

X 

0 

X 

X 

0 

X 

0 

0 

X 

X 

0 

X 

X 

0 

X 

0 

X 

0 

0 

X 

X 

0 

X 

X 

0 

0 

X 

0 

X 

X 

0 

X 

X 


(108) 


(109) 


( 110 ) 


( 111 ) 


( 112 ) 


(113) 


Since Mr is diagonal in this case, the constraints mainly come from four zeros in 
and they are consistent with the results discussed in Refs. [7, 28]. 


36 



A.3 Nr = 2 

If Mji has two zeros, there are four structures of Mji that are not equivalent under permu¬ 
tation. They are 


Mrr — 


Mrc — 


0 

0 

X 


0 

X 

X 

0 

X 

X 

, Mri, — 

X 

X 

0 

X 

X 

X 


X 

0 

X 

X 

0 

0 


X 

X 

X 

0 

X 

X 

, Mr^ = 

X 

0 

X 

0 

X 

X 


X 

X 

0 


(114) 


All other structures of Mr can be obtained by a permutation of the columns and rows of 
the above structure. If Mr, have more than 6 zeros, M or C would have at least three zeros 
for all three structures of Mr, so here we only consider Nr, < 6. 


A.3. 

1 . 


1 Mjia 

Nr) = 6. Of the 84 different structures of Mr in this case, 78 lead to at least three 
zeros in M or C, and 6 lead to Class X constraints for M. 


2. Nr, = 5. Of the 126 different structures of Mr, in this case, 78 lead to at least three 
zeros in M or C, 18 lead to Class X constraints on M, 6 lead to one diagonal zero in 
M, and 6 lead to one diagonal zero in C. In addition, there are 6 structures that lead 
to Class 2 nonstandard constraints: 



0 

0 

X 


0 

0 

X 

Mr = 

X 

0 

X 


X 

0 

0 


X 

0 

0 


X 

0 

X 


which lead to the constraints, Mge = 0 and det M = 0; 



X 

0 

X 


X 

0 

0 

Mr = 

0 

0 

X 


0 

0 

X 


X 

0 

0 


X 

0 

X 


(115) 


(116) 
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which lead to the constraints, = 0 and det M = 0; 



X 

0 

X 


X 

0 

0 

Md = 

X 

0 

0 


X 

0 

X 


0 

0 

X 


0 

0 

X 


which lead to the constraints, M^-j. = 0 and det M = 0. 

Also, there are 6 structnres lead to Class 6 nonstandard constraints: 



0 

X 

0 


0 

X 

0 

Md = 

X 

0 

X 


X 

0 

0 


X 

0 

0 


X 

0 

X 


which lead to the constraint, M^eCee = detM; 



X 

0 

X 


X 

0 

0 

Md = 

0 

X 

0 


0 

X 

0 


X 

0 

0 


X 

0 

X 


which lead to the constraint, = det M; 



X 

0 

0 


X 

0 

X 

Md = 

X 

0 

X 


X 

0 

0 


0 

X 

0 


0 

X 

0 


which lead to the constraint, MrrCrr = det M. 

The remaining 6 structnres lead to the Class 9 constraint, det M = 0: 


X 

X 

0 


X 

X 

0 


X 

0 

0 

X 

0 

0 


0 

X 

0 


X 

X 

0 

0 

X 

0 


X 

0 

0 


0 

X 

0 

X 

0 

0 


0 

X 

0 


0 

X 

0 

0 

X 

0 


X 

X 

0 


X 

0 

0 

X 

X 

0 


X 

0 

0 


X 

X 

0 


( 117 ) 


(118) 


(119) 


( 120 ) 


( 121 ) 


( 122 ) 
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A.3.2 Mm 

1. Nd = 6 . Of the 84 different structures of Mu in this case, 78 lead to at least three 
zeros in M or C, and 6 lead to Class Y constraints on C. 


2. No = 5. Of the 126 different structures of Mo in this case, 78 lead to at least three 
zeros in M or C, 24 lead to one diagonal zero in C, and 12 lead to one off-diagonal 
zero in C. The remaining 12 structures lead to the Class 9 constraint, det M = 0: 


X 

X 

0 


X 

X 

0 


X 

0 

0 

X 

0 

0 


0 

X 

0 


X 

X 

0 

0 

X 

0 


X 

0 

0 


0 

X 

0 

X 

0 

0 


0 

X 

0 


0 

X 

0 

0 

X 

0 


X 

X 

0 


X 

0 

0 

X 

X 

0 


X 

0 

0 


X 

X 

0 

X 

0 

X 


X 

0 

X 


X 

0 

0 

X 

0 

0 


0 

0 

X 


X 

0 

X 

0 

0 

X 


X 

0 

0 


0 

0 

X 

X 

0 

0 


0 

0 

X 


0 

0 

X 

0 

0 

X 


X 

0 

X 


X 

0 

0 

X 

0 

X 


X 

0 

0 


X 

0 

X 


(123) 


(124) 


(125) 


(126) 


A.3.3 Ma, 

1. No = 6 . Of the 84 different structures of Mo in this case, 78 lead to at least three 
zeros in M or C, and 6 lead to a block diagonal matrix for M. 


2. No = 5. Of the 126 different structures of Mo in this case, 72 lead to at least three 
zeros in M or C, 12 lead to a block diagonal matrix for M, 24 lead to one off-diagonal 
zero in M, and 12 lead to one off-diagonal zero in C. The remaining 6 structures lead 
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to the Class 9 constraint, det M = 0: 


0 

X 

X 


0 

X 

X 


0 

X 

0 

0 

X 

0 


0 

0 

X 


0 

X 

X 

0 

0 

X 


0 

X 

0 


0 

0 

X 

0 

X 

0 


0 

0 

X 


0 

0 

X 

0 

0 

X 


0 

X 

X 


0 

X 

0 

0 

X 

X 


0 

X 

0 


0 

X 

X 


( 127 ) 


(128) 


A.3.4 Mud 


1. Nd = 6. Of the 84 different strnctures of Mr, in this case, 78 lead to at least three 
zeros in M or C, and 6 lead to Class Z constraints on M. 


2. Nd = 5. Of the 126 different structnres of Mr, in this case, 84 lead to at least three 
zeros in M or C, 12 lead to Class Z constraints on C, and 24 lead to one diagonal zero 
in C. The remaining 6 strnctnres lead to the Class 9 constraint, detM = 0: 


0 

X 

X 


0 

X 

X 


0 

X 

0 

0 

X 

0 


0 

0 

X 


0 

X 

X 

0 

0 

X 


0 

X 

0 


0 

0 

X 

0 

X 

0 


0 

0 

X 


0 

0 

X 

0 

0 

X 


0 

X 

X 


0 

X 

0 

0 

X 

X 


0 

X 

0 


0 

X 

X 


(129) 


(130) 


A.4 Nji — 1 


For Nd = 6, of the 84 different strnctnres of Mr, in this case, 78 lead to more than three 
zeros in M or C, and 6 lead to one diagonal or off-diagonal zero in C. 
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